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We review selected results for Higgs boson production at Linear Colliders in the framework of the general 
^< , Two-Higgs-Doublet Model (2HDM). We concentrate on the analysis of i) the pairwise production of 

Q\ ' neutral Higgs bosons (h°A ,H°A ); and ii) the neutral Higgs boson- strahlung modes (h°Z°, H°Z°). We 

identify sizable production rates, in the range of a ~ O(10 — 100) fb for *J~s = 0.5 TeV, alongside 
with large quantum effects (5 r ~ ±50%), which we can fundamentally track down to the enhancement 
power of the triple-Higgs self-interactions. This constitutes a telltale signature of the 2HDM, with no 
counterpart in e.g. the Minimal Supersymmetric Standard Model (MSSM). We compare these results 
£1h with several complementary double and triple 0(al w , ai w ) Higgs-boson production mechanisms and 

spotlight a characteristic phenomenological profile which could eventually be highly distinctive of a non- 
supersymmetric two-Higgs-doublet structure. 
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> ' 1 Introduction 

Deciphering the fundamental nature of Electroweak Symmetry Breaking (EWSB) lies at the very frontline of 
the current theoretical and experimental research in Particle Physics. Even in spite of the tantalizing Higgs 
boson candidates recently identified by the ATLAS and CMS experiments pQ, a long way might yet stand 
ahead of us until we are able to convincingly close in on such a longstanding conundrum. In particular, were 
^vq . these signatures finally confirmed, a first question to be answered would be whether they can be described 
within the minimal framework of the Standard Model (SM) or, on the contrary, if they should rather be 
attributed to an extended EWSB sector [2J. One canonical example of the latter is the general Two-Higgs- 
Doublet Model (2HDM) [3]. The model is built upon a second scalar SUl(2) doublet with Y = +1 weak 
hypercharge. This results into a larger Higgs boson spectrum of five physical Higgs fields: neutral CP-even 
(h°,H°), neutral CP-odd (A ) and charged H ± . Such a simple, but yet non-minimal extension of the SM 
Higgs sector has gathered growing attention over the past years [3] and become a cherised setup for model 
builders and phenomcnologists alike. Besides the many novel, and usually highly distinctive features put 
forward by the model in multifarious domains - from collider to flavor physics or astrophysics -, the 2HDM 
provides a suitable low-energy description to many UV completions of the the EWSB dynamics. 

The model is fully specified once we fix i) the masses of the physical Higgs bosons; ii) the ratio tan j3 = 
(H®) / (Hf) of the two Vacuum Expectation Values (VEVs) giving masses to the up- and down-like quarks; 
iii) the mixing angle a between the two CP-even states; and iv) one remaining Higgs boson sclf-couplingin 
the Higgs potential, hereafter dubbed A5. We note in passing that the Higgs sector of the MSSM [5] 
corresponds to a particular realization of the general (unconstrained) 2HDM, for which the invariancc under 
SUSY transformations enforces a number of additional restrictions - most significantly, the Higgs boson 
self-interactions become tied to the gauge couplings. This situation is remarkably different in the general 
2HDM, where the size of these self-interactions is fundamentally unrestricted and it only becomes limited, 
in practice, by the interplay of theoretical consistency conditions (unitarity [BJ, vacuum stability [7]) and 
experimental bounds (viz. the excluded Higgs boson mass ranges from the from direct collider searches, 
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Figure 1: Sample of one- loop Higgs- mediated Feynman diagrams which account for the bulk of the quantum 
effects to the neutral Higgs pair e + e~ — > h°A° (upper row) and the Higgs-strahlung mechanisms c + c~ — > h°Z° 
(lower row). 



and also the constraints derived from electroweak [SJ and flavor physics observables [5]). A detailed account 
of these restrictions and of the model setup can be found in Ref. [TO]- For comprehensive analyses of the 
2HDM parameter space constraints, see e.g. Refs. [TTj . 

Following the eventual discovery of the Higgs boson(s) at the LHC, it will be crucial to address the precise 
experimental determination of the corresponding quantum numbers, mass spectrum and couplings to other 
particles. A linear collider (linac) can play a central role in this enterprise [T^]. Dedicated studies have 
exhaustively scrutinized the phenomcnological imprints of the basic 2HDM Higgs boson production modes, 
such as e.g. i) triple Higgs, e + e~ — > 3h [13]; ii) inclusive Higgs-pair through EW gauge boson fusion, e + c~ —> 
V*V* —> 2h + X [14]; iii) exclusive Higgs-pair e + e~ — > 2h [TOJUH]; an d i y ) Higgs strahlung, or associated 
production with a weak gauge boson e+c~ hV [Tg], with \h = h°,A°,H°,H ± ] and [V = Z°,W ± jU. As 
a common highlight, all these studies report on sizable production rates and large quantum effects, arising 
from the potentially enhanced Higgs self-interactions. Interestingly enough, Higgs boson searches at e + e~ 
colliders may also benefit from alternative running modes, particularly from 77 scattering. Processes such 
as 77-induced production of single (77 — > h) and double (77 — > 2/i) Higgs bosons have been studied from 
this viewpoint. These entirely operate at the quantum level, via an effective (loop-induced) Higgs/photon 
interaction <? 77 /j that we may regard as a direct probe of non-standard (charged) degrees of freedom coupled 
to the Higgs sector. The aforementioned single Higgs channels have been considered in the framework 
of the SM [IS], the 2HDM [TH] and the MSSM gUHH] and are known to exhibit excellent experimental 
prospects, not only due to the clean environment inherent to a linac machine, but also owing to the high 
attainable 77 luminosity, and the possibility to tune the 7-beam polarization as a strategy to enlarge the 
signal- versus-background ratios^. 

2 Numerical analysis 

In this contribution we review two particular 2HDM Higgs boson production modes at a linear collider, to 
wit: i) the pairwise production of neutral Higgs bosons e + e~ — > h A°/H°A ; and ii) the associated production 
of a neutral Higgs and a Z° bosons, e + e _ — > h Z°/H°Z - the so-called Higgs-strahlung mechanism, which 
we can regard as the 2HDM analog(s) to the Bjorken process in the SM The motivation herewith is 

threefold: i) a first focus point is to seek for the most favorable regions across the 2HDM parameter space, 
for which the Higgs boson production rates become optimal, and to correlate them to alternative multi-Higgs 
production modes; ii) second, we aim at quantifying the importance of the radiative corrections associated 
to these processes; iii) and third, we shall examine the impact of the Higgs boson self interactions and their 

1 For related work in the context of MSSM Higgs boson production see e.g. I17| . 
2 Analoguc studies for the 77 — > hh mode are available e.g. in Ref. 1221 . 
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Figure 2: Tree-level and loop-corrected cross section [cr(e + e~ — > h°A )] (in fb), alongside with the relative 
size of the radiative corrections [5 r = [<r' 0+1 ' — er(°']/er ] (in %), as a function of y^s (left, center) and A5 
(right). We fix tan/3 — 1.2 (compatible with the lower tan/3 bound from B® — B® data j9]) and examine 
the representative choices a = f3 (maximum gh°A°z° tree- level coupling) and a = tt/2 (fcrmiophobic limit 
for h° in type-I 2HDM). The influence of the Higgs self-interactions is assessed by dialing the value of the 
parameter A5. The shaded areas on the left (rcsp. right) are excluded by unitarity (resp. vacuum stability). 



potentially enhanced strenght. The leading order production rates merely depend on the Higgs couplings to 
the Z boson. In other words, they are entirely subdued by the gauge symmetry - and hence they feature no 
disclosing scenarios between the general 2HDM and e.g. the MSSM. The resulting phenomenological portray, 
however, may clearly depart once the quantum effects to such couplings are considered. Vertex corrections, 
in particular, turn out to be sensitive to the triple Higgs self-interactions through the interchange of virtual 
Higgs bosons which are then linked to the external Higgs boson legs. A sample of such Higgs-mediated 
one-loop diagrams is displayed in Fig. [T] These effects we can roughly describe by a loop-induced form 
factor, which spells out how the strength of the bare Higgs-to-gauge boson couplings is modified: 



5/iA°Z° 9hA°Z° 



HHH I 



16tt 2 



f{Ml/s,Ml /s) 



(1) 



Here Xhhh stands for generic triple Higgs self-interaction, and 1/16-7T 2 for the standard loop integral sup- 
pression factor. By /(M 2 /s, M? /s) we denote a generic rational function involving the ratios of the different 
mass scales taking part in the process. The above expression (p} indicates how the Higgs-to-gauge boson cou- 
plings, which are entirely anchored by the gauge symmetry at the leading-order, may be strongly promoted 
at one-loop through the indirect effect of the Higgs boson self-couplings - unlike their MSSM counterparts. 



2.1 Calculational setup 

Throughout our study we make use of the standard algebraic and numerical packages FeynArts, Form- 
Calc and LoopTools Updated experimental constraints (from EW precision data, low-energy flavor- 
physics and the Higgs mass regions ruled out by direct collider searches), as well as theoretical consistency 
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Figure 3: Contour plots for the relative size of the one-loop quantum corrections 5 r = [er( n+1 ' — a^j/a (in 
%) to the e + e~ — ^ h°A° (left panel) and e + e~ — > h°Z° (right panel) total rates, as a function of tan /3 and A5. 
We fix a = /3 (for e + e~ — > h°A°) and a = /3 — n/2 (for c + c~ — > h°Z°), in which cases their respective born- 
level couplings maximize. For the Higgs boson masses we use Set A from Tab. [T] The linac center-of-mass 
energy is taken to be >/s = 500 GeV. The shaded areas in the upper (resp. lower) patches of the tan/3 — A 5 
plane are excluded by unitarity [5] (resp. vacuum stability [7]) bounds. The vertical grey strip accounts for 
the lower limit tan/3 ~ 1.18 stemming from — data [S]- 



conditions (perturbativity, unitarity and vacuum stability) are duly taken into account [51 171 151 [TT1 1251 
For definiteness, we set along two Higgs boson mass benchmark choices A and B, as quoted in Tab.[TJ 

M h o [GeV] M H o [GeV] M A o [GeV] M H ± [GeV] 



Set A 
Set B 



130 
130 



200 
150 



260 
200 



300 
160 



Table 1: Choices of Higgs boson masses employed throughout our calculation. 



In order to carry out the complete one-loop computation we are entitled to define suitable UV coun- 
terterms, in particular for the renormalization of the Higgs boson masses and wave functions. These we 
can express in terms of the Higgs 2-point functions at order 0(a ew ). Conventional on-shell renormalization 
conditions - see e.g. Ref. [27] - are extended to the 2HDM Higgs sector. In particular, the relations 



ReS^o A o(<? 2 ) 



q 2 = M' 



0; ReE A o z o(g 2 ) 



AO 



9 2 = A'Io 



0, 



anchor the wave function renormalization of the Higgs doublets, and thereby of all the physical Higgs fields. 
The remaining Higgs boson masses, as well as the mixing angle a, are determined via on-shell conditions 
imposed on their respective self-energies (including the h° — H° kinetic mixing). The parameter tan/3 is fixed 
via Eq. alonside with one additional condition on the Higgs boson tadpoles. An exhaustive account of 
the renormalization procedure is available in Ref. |10j . 
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Figure 4: Tree-level and loop-corrected cross section [<r(e + e — s- h°A )] (in fb) and relative size of the 
radiative corrections [S r ] (in %), as a function of y/s (left, center) and A5 (right), for an equivalent setup to 
that of Fig. [2J 



2.2 Higgs boson pair production at 0(pP ew ): e+e — > h A°/H°A 

For definiteness, we focus on the light Higgs mode [h°A°] and specialize our results for the Higgs mass 
spectrum defined by Set A of Tab. [TJ We quantify our analysis by means of i) the Born-level, [cr*- - 1 ], and 
1-loop cross sections, [ct' 0+1 ^] - in which we include the full set of 0{a\ w ) corrections, and also the leading 
0(a^ w ) ones which arise from the squared of the scattering amplitude [10]; ii) the relative size of the 1-loop 
corrections, via the parameter S r = [o"( 0+1 ) — a^]/a^°\ The upshot of our findings, as summarized in Fig.[2J 
highlights substantial production rates, which fall roughly in the range of 2 — 15 fb for y/s = 0.5 TeV - this 
is, up to barely 10 3 — 10 4 events per 500 fb -1 ; and eventually very large quantum corrections, of the order 
\Scr\/cr ~ 20 — 60%, which can be either positive (for s/s ~ 0.5 TeV) or negative ( s/s > 1 TeV) and fairly 
independent on the details of the Higgs mass spectrum, the particular value of the tree-level coupling [<7^,a°z] 
and the actual type of 2HDM under consideration - namely, whether we specifically target at the type-I 
or II realizations of the 2HDM. The evolution of a^ 0+1 ^ and S r as a function of y/s for different A5 values 
evinces how critically the quantum effects depend on the Higgs self-interaction enhancements. The quadratic 
dependence on the parameter A5, a ~ (a — 6A5) 2 , as shown in the rightmost panel of Fig.[5J nicely illustrates 
the dominance of the Higgs mediated one-loop diagrams - these are indeed sensitive to the product of two 
triple Higgs self-interactions. As a complementary viewpoint, in the left panel of Fig.[3]we display the profile 
of the radiative corrections 5 r to the total cross-section [(T(h°A )] along the tan j3 — A5 plane, again for Set 
A of Higgs boson masses, a = /3 and a linac center-of-mass energy of y/s = 0.5 TeV. The choice a = (3 
maximizes the tree-level cross section. Unitarity and vacuum stability limits [7] (lower and upper shaded 
areas, respectively) restrict the largest attainable quantum effects to regions with tan/3 ~ 1 and | A5 1 ~ 5 — 10 
(A5 < 0). The central grey band depicts the lower limit tan/3 ~ 1.18 ensuing from — B% [9]- 

2.3 Associated Higgs/Z°-boson production at 0(a 3 ew ): e+e" ->■ h Z°/H°Z 

Again, without loss of generality, we concentrate on the light Higgs mode [h°Z°] and arrange the mass 
spectrum as in Set A of Tab. [1] Our results are shown in FigH] In this case we obtain typical cross sections 



5 



Process 


cr(V5=0.5TeV)[fb] 


cr(Vs = 1.0TcV)[fb] 


a(y/s = 1.5TeV)[fb] 


h (, A J 


26.71 (5 r = 31.32%) 


4.07 


1.27 


H°Z° 


19.09 (S T = -61.56%) 


3.73 


1.47 


h°H°A° 


0.02 


5.03 


3.55 


H°H+H~ 


0.17 


11.93 


8.39 


h°h° + X 


1.47 


17.36 


38.01 



Table 2: Compared cross section (in fb) for different associated, pairwise and triple Higgs boson production 
mechanisms at 0(a^ w , a^ w ), for tan/3 = 1, a = (3 and A5 = —10. The Higgs boson mass spectrum we fix as 
in Set B of Tab. [TJ The relative size of the one-loop corrections [S r ] for the Higgs pair and Higgs strahlung 
mechanisms is quoted in brackets. 



in the range of a(h°Z°) ~ O(10 — 100) fb, with very significant (and systematically negative) radiative 
corrections (up to order S r ~ —50%), reaching their maximum again in the parameter space regions with 
tan/3 ~ 0(1) and | A5 1 ~ O(10). Such a characteristic pattern of negative quantum effects we can relate 
to the dominance of the finite wave function corrections to the external Higgs boson fields - this being the 
only contribution which retains a quadratic dependence on Xhhh at one loop, as we can also read off the 
rightmost panel of Fig. 21 The relative size of the quantum effects [8 r ] and its interplay with the relevant 
constraints is examined in the right panel of Fig. [3] as we move across the tan j3 — A5 plane. Set A of Higgs 
boson masses, a fixed value of a = (3 — ir/2 and a linac center-of-mass energy to y/s = 0.5 TeV are employed 
throughout. Worth noticing is that the S r isocurves are not responsive to a change of tan/3. This follows 
directly from the analytical structure of all the relevant couplings for the particular setup a = (3 — it/ 2 [TB] 
- which corresponds to the decoupling (or SM-like) limit of the 2HDM. 

3 Discussion and closing remarks 

Higgs boson self-interactions constitute a paradigmatic structure of extended Higgs sectors of non-supcrsymmetric 
nature. The general (unconstrained) Two-Higgs-Doublet Model is a canonical example of the latter. Here, 
the triple and quartic Higgs boson self-interactions are not subdued by the gauge symmetry. This entails 
two major consequences, which are in stark contrast to analogue models, such as e.g. the MSSM: i) the 
Higgs boson spectrum is fully unconstrained; this is to say, no limitations on the mass splittings between the 
physical Higgs boson fields must be assumed a priori; ii) by the very same token, the Higgs boson self interac- 
tions are also fundamentally unrestricted, and hence may accomodate sizable enhancements. Both features 
are tamed in part by stringent theoretical and phenomenological constraints (unitarity, vacuum stability, 
electroweak precision and flavor physics) but nevertheless open up a plethora of rich, and highly distinc- 
tive, phenomenological possibilities. So much so, the analysis of collider observables which are sensitive to 
the Higgs self-interactions, cither directly or through quantum corrections, may bring forward instrumental 
handles to disclose non-SUSY vs SUSY multi-doublet Higgs structures. 

In this contribution we have concisely revisited two particular examples of Higgs boson production 
from e + e~ colliders within the 2HDM context, these are the pairwise production of neutral Higgs bosons 
(h A°/H°A ) and the Higgs-strahlung channels (h Z°/H°Z ). We have portrayed their phenomenology at a 
future linac and have spelled out the features that singularize the 2HDM scenarios with large Higgs boson 
self-couplings. Our findings can be outlined as follows: 

• Large Higgs boson production rates, in the ballpark of a2h,hz ~ O(10 — 100) fb for y/s = 500 

GcV, and yet of few dozens of fb in the TcV-range center of mass energies - this would correspond to 
rates of O(10 2 — 10 5 ) events for an integrated luminosity of 500 fb -1 . 

• Large quantum effects, which may reach up to S r = [o"( 0+1 ' — a^]/cr^ ~ ±50%, preferably 
realized within the tan/3 ~ 0(1) and |Ag| ~ O(10) domains of the 2HDM parameter space. These may 
alternatively lead to characteristic enhancements (e.g. for a(2h) at y/s ~ 500 GeV), or suppressions 
(for <r(hZ) and also for a(2h) at larger y/s). 
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• A generic phenomenological pattern, in the sense that the above observations barely depend on 
the very choice of Higgs masses nor the type of Yukawa couplings to fermions, and they hold for broad 
regions across the tan j3 — sin a plane. 

Interestingly enough, enhancements of the Higgs boson production rates have also been put forward in the 
literature for alternative multi-Higgs production processes, such as the e + c~ — > hhh [13] and e + e~ — > VV* —> 
hh + X [14] channels. In this vein, let us consider the following choice of parameters: tan/3 = 1, a = (3 and 
A5 = —10, along with Set B of Higgs boson masses from Tab.Q] This particular configuration saturates the 
unitarity bounds, and thus maximizes the impact of the Higgs boson self-interactions. If we now combine 
the evaluation of the total production rates for all those different production channels, we conic up with the 
cross-correlated set of predictions displayed in Tab. [5] These results reflect the great complementary of the 
different multi-Higgs states at different center-of-mass energies. Likewise, the correlation of large negative 
quantum effects on the Higgs-strahlung channels with the presence of significant positive (for yj~s < 500 
GcV) and negative (for y/s > 600 GeV) quantum effects on the double Higgs production may eventually 
constitute an additional hint at a generic (unconstrained) 2HDM dynamics. Notice once more that, in all 
these cases, the reported pattern of signatures crucially relies on the strenght of the 3H self-couplings. No 
analogue picture could then be attributed e.g. to the MSSM. 

On balance, there is little doubt that a linear collider qualifies as a most cherised tool to carry to 
completion the Higgs boson research program. Owing to its superbly clean environment, a linac facility 
should enable accurate measurements of the Higgs boson masses, gauge and Yukawa couplings, as well as of 
the Higgs boson self-interactions themselves. This means, it could provide us with the firmest possible grip 
on the fundamental structure of the EWSB sector. In this context, our results underline the possibilities of 
the Higgs boson self-interactions as a trademark dynamical feature of the generic 2HDM. We prove their 
capabilities to rubber-stamp significant - and highly distinctive - fingerprints on multi-Higgs production 
processes, either at the leading order or through quantum corrections, and conclude that these might well 
constitute a pristine window towards non-standard, non-supcrsymmctric Higgs sectors. 
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